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The effect of continuum couplings in the fusion of the halo nucleus ^^Be 
on ^''^Pb around the Coulomb barrier is studied using a three-body model 
within a coupled discretised continuum channels (CDCC) formalism. We in- 
vestigate in particular the role of continuum-continuum couplings. These are 



O ' found to hinder total, complete and incomplete fusion processes. Couplings to 

^ ■ the projectile lpi/2 bound excited state redistribute the complete and incom- 

^ , plete fusion cross sections, but the total fusion cross section remains nearly 

constant. Results show that continuum-continuum couplings enhance the 



irreversibility of breakup and reduce the flux that penetrates the Coulomb 
barrier. Converged total fusion cross sections agree with the experimental 
ones for energies around the Coulomb barrier, but underestimate those for 
energies well above the Coulomb barrier. 



Introduction: The existence and the role of the breakup process of weakly bound pro- 
jectiles in complete fusion and scattering mechanisms have been extensively investigated 
O ' in recent years both theoretically and experimentally but there is not yet any 

^ . definitive conclusion. There are contradictory theoretical works which predict either the 
suppression or the enhancement |^ of the complete fusion cross section due to the 

CJ I coupling of the relative motion of the nuclei to the breakup channel. 

^ ■ Recent coupled channels calculations for ^^Be-l-^^'^Pb [|| have shown that the coupling 

of the relative motion to the breakup channel has two effects, depending on the value of the 
bombarding energy, namely (i) a reduction of the complete fusion cross sections at energies 
^ i above the Coulomb barrier due to the loss of incident flux, and (ii) an enhancement of the 
I complete fusion cross sections at energies below the Coulomb barrier due to the dynamical 
renormalisation of the nucleus-nucleus potential. Using the iso centrifugal approximation 
and an incoming boundary condition inside the barrier, this calculation did not include the 
effect of the projectile's halo structure on the monopole projectile-target potential. Nor did 
it include the excitation to partial waves other than in the continuum, or the continuum- 
continuum and bound excited states couplings in either reaction partner. Moreover, only a 
small interval of energy for continuum states (up to 2 MeV) was considered. 

The couplings between continuum states have been shown to be crucial to understand 
the breakup of ®B on a ^^Ni target at low energy Eiab = 25.8 MeV JT^pD| . Therefore, 



it could be expected that continuum- continuum couplings significantly affect the role of 
breakup process in fusion of halo nuclei around the Coulomb barrier. We believe that 
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continuum-continuum couplings enhance the irreversibihty of the breakup process (thus 
once the projectile gets dissociated, it will find it very hard to find its way back to the 
bound states). Thus, the continuum-continuum couplings are expected to reduce the fiux in 
bound projectile channels and, therefore, should inhibit (at least) the complete fusion. We 
expect that, with an increasing breakup subspace, governed by both the maximum energy 
and continuum partial waves, continuum-continuum couplings will reduce the fusion cross 
sections until convergence is reached. 

The aim of this paper is to clarify the role of these continuum couplings on the fusion 
of the halo projectile ^^Be on a ^o^Pb target at energies around the Coulomb barrier. Cal- 
culations are carried out using a three-body model [l^JSn] in the framework of the CDCC 



formalism |l21| , |22| . Full coupled channels calculations are performed with the code FRESCO 



In calculation of fusion cross sections, we simultaneously include (i) the effect of the 
projectile's halo structure on the projectile-target potential, (ii) both the transition to its 
bound excited state and its dissociation caused by inelastic excitations to different partial 
waves in the continuum, induced by the projectile fragments-target interactions (coulomb 
-|- nuclear), and (iii) couplings (bound-continuum and continuum-continuum) between its 
excited states. We do not include transfer or inelastic channels of the target. Fusion cross 
sections for projectile bound channels and for projectile breakup channels will be defined in 
terms of a short-ranged imaginary bare potential defined in the center of mass coordinate of 
the projectile in conjunction with channel dependent wave functions for the projectile-target 
radial motion. 

Model: In the ^^Be + ^'^^Pb reaction, the three-bodies involved are the ^°Be core (C), the 
valence (halo) neutron (v), and the ^osp]^ target (T). 

Let R be the coordinate from the target to the center of mass of the projectile, and "r^ 
the internal coordinate of the projectile. The position coordinates of the projectile fragments 
with respect to the target are T^vt = R + and Rct = R — ■, where is the 

mass of the projectile. 

The dynamics of the three-bodies is described by the Schrodinger equation in the over-all 
centre of mass system. The two-body potentials Vct, Kt and depend only on relative 
coordinates indicated as their arguments and do not excite the internal degrees of freedom of 
the core and the target nucleus. Following 0, these potentials are considered as real, but in 
addition we use for fusion a bare short-ranged (well inside the Coulomb barrier) imaginary 
central potential iWF{R) defined in the center of mass of the projectile for the projectile- 
target radial motion. The use of this short-ranged imaginary potential is equivalent to the 
use of an incoming boundary condition inside the Coulomb barrier to study fusion P^JT^ . 



In order to describe the breakup of a projectile such a ^^Be, we consider the inelastic 
excitations induced by Vct-, Kt in the n + ^°Be system from the ground state (l>(is)j ni''^) 
excited states in the continuum U(is)j^k{f), for some wave- number k and partial wave /, and 
also couplings between all such continuum states. The use of such single energy eigenstates, 
however, would result in calculations of form factors for continuum-continuum couplings 
which do not converge, as the continuum wave functions are not square integrable. The 
CDCC method [PT| , P^ is used to obtain square integrable continuum bins states <P{is)j,[ki,k2](^) 
averaged over a narrow range of wave- numbers [ki,k2]- We label these bin states by their 
wave-number limits [ki, /C2] and their angular momentum quantum numbers {ls)j. The 
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bound states of the projectile (p{is)j,nir) and the single energy scattering wave functions 
U(is)j,kij) which form the continuum bins 0(is)j,[fci,fc2](^)) are obtained by solving a Schrodinger 
equation with the potential V^(j which may be /-dependent. The bin wave functions are 
defined as 



(Zs)i,[fci,fc2]l^J 



TlN 



k2 



w{k)e '^''U(is)j,k{r)dk, 



(1) 



where 5k is the scattering phase shift for 'U(;s)j,fc(r). The normalisation constant is = 

k2 

J \w{k)\'^dk for the assumed weight function w{k), here taken to be either unity for non-s- 

ki 

wave bins or k for s-wave bins. These bin states are normalised {4>\4>) = 1 once a sufficiently 
large maximum radius r;,j„ for r is taken. They are orthogonal to any bound states, and are 
orthogonal to other bin states if their energy ranges do not overlap. The phase factor e"*"^* 
ensures that they are real valued for real potentials V^^^. 

The radial wave functions faj{R) for the projectile-target relative motion satisfy the set 



of coupled equations |T9 | 




L{L + 1] 

i?2 



+ VlJR)+iWF{R) + e^-E 



faAR) 



= E ^''-'vl,AR)fo.'j{R), (2) 

where /x, L, J, E and a {{L,l, s, j,n or [A;i,/c2]}) denote the projectile-target reduced mass, 
the projectile orbital angular momentum, the total angular momentum, the total energy, 
and the set of quantum numbers, respectively. For unbound states of the projectile, ea is 
the mean energy of continuum bin [ki, ^2] , or ea < for bound states. V^.^, describes the 
coupling between the different internal states 0q-(^) of the projectile 

yL'0) =< M^) I VctC^ct) + V,t{^,t) I 0a'( r^) > . (3) 

Assuming that the potentials Vct and V^t are central, the Legendre multipole potentials 
can be formed as 

1 y - 

AK{R,r) = -J [Vct{Rct) + V,T{^,T)]PK{x)dx, (4) 
-1 

where K is the multipole and x = r ■ R is the cosine of the angle between "P" and R . Since 
the spin s of the neutron is fixed, the coupling form factor (|^) between states 0Q,'(r) and 
(pair) is 

VL'iR) = Y.i-^y^''''"T3'U'LL'(2K + l)W{33'll'- Ks)W{jj'LL'; KJ) 

K 

''(foo)(foo) jMr)AK{R,r)(Pa'{r)dr. (5) 
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Eqs. (g) are solved with the usual scattering boundary conditions . 

The total fusion cross section atot is defined in terms of that amount of fiux which leaves 
the coupled channels set because of the short-ranged imaginary potential iWf^R)- 

Since complete fusion is a process where all the nucleons of the projectile are captured 
by the target nucleus and following 0, we define in our model the complete fusion cross 
section as the absorption cross section from projectile bound channels (complete fusion from 
both ground state (elastic) and bound excited state) 

<'cF = ^Y.i2J+l)Pj, (6) 

where E is the bombarding energy and Pj is the complete fusion probability for the partial 
wave J. The complete fusion probability Pj is 



/i(2E//i)i/2 



oo 

E / I fo^AR) r {-WF{R))dR. (7) 



a {ea< 0) 



The complete fusion cross section (|^)-(0) represents a lower limit for the physical complete 
fusion cross section, since we have assumed no capture of all projectile fragments (^°Be and 
the halo neutron) from breakup channels. In reality, these events should contribute to the 
complete fusion, but cannot be distinguished in our model from the capture of only one 
projectile fragment. 

The incomplete fusion aicF (fusion of ^°Be) is then defined as the absorption from 
breakup channels 

Results and discussion: The experimental spectrum of ^^Be exhibits a 1/2"'" ground state 
and a single, l/2~, bound excited state with energies of -0.50 MeV and -0.18 MeV, respec- 
tively. In a pure single-particle picture, the ground and the bound excited states of ^^Be 
have 2si/2 and lpi/2 single-particle configurations, respectively. These configurations can 
be associated with single-particle states generated by different V},(j Woods-Saxon potentials 
p5| including a spin-orbit term. For the 2si/2 state, we use a Woods-Saxon potential with 
parameters Vq = -51.51 MeV, = 1.39 fm and a = 0.52 fm. For the lpi/2 state, we use 



a Woods-Saxon potential including a spin-orbit term, similar to that used in p5|, with the 
same geometry, i.e. Vq = -30 MeV, tq = 1.39 fm, a = 0.52 fm and V^*'°' = 4.39 MeV. 

First, we study qualitatively the effect of continuum couplings on fusion cross sections 
by using a reduced breakup subspace with regard to the maximum energy of the projectile 
continuum states. A continuum breakup subspace with partial waves Si/2, P1/2, P3/2, d3/2 
and d5/2, for the halo neutron-^°Be core relative motion, is used. For each partial wave, the 
continuum subspace is discretised in 6 bins which are equally spaced in wave-number k, up 
to a maximum wave-number kmax = 0.3612 fm~^ (a maximum energy of 3 MeV), with a step 
of Ak = 0.0602 fm~^. In Fig.l, we illustrate the continuum discretisation used to define the 
energy bins included in these calculations. The calculation is thus performed with 30 excited 
continuum channels. The s- and p-wave continuum states have been consistently generated 
by the same potential V^q as that of the bound state of the same angular momentum /. The 
d-wave continuum states have been generated by the same potential as that of the p-waves. 
In the present work, Woods-Saxon parametrisations given in |]2^ and in are used for 



the nuclear part of the potentials Vct (Vq = -46.764 MeV, tq = 1.192 fm and a = 0.63 fm) 
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and VvT (Vq = -44.019 MeV, tq = 1.27 fm and a = 0.67 fm), respectively. A short-ranged 
Woods-Saxon potential Wp with parameters Vq = -50 MeV, tq = 1 fm and a = 0.1 fm is used 
for the fusion potential. The results depend only weakly on the geometry of this potential, 
as long as it is well inside the Coulomb barrier and strong enough that the mean-free path 
of the projectile inside the barrier is much smaller than the dimensions of Wp. The fusion 
cross sections for Vq = -50 MeV are those for Vq = -10 MeV changed by ~ 1%. 

Since we are interested in fusion cross sections, partial waves for the projectile-target 
relative motion up to only L^ax = 50 (partial-wave total fusion cross section ~ 10"'^ mb) are 
included. Our calculations include monopole, dipole and quadrupole contributions {K=0,1 
and 2) of the potentials Vct and V^t for both nuclear and Coulomb parts. The couplings 
V^.^,{R) are taken into account up to a projectile-target radial distance Rcoup = 100 fm. To 
calculate both the continuum bins (|1|) and couplings V^.^,{R) (^, which include these bins, 
radii r < rf,j„ = 100 fm are used. 

Figs. 2a and 2b show fusion cross sections as a function of the bombarding energy in 
the center of mass system. For comparison, we present cross sections in the absence of 
couplings (thin solid curve). In Fig. 2a, calculations include transitions from and to the 
projectile bound states, but do not include continuum-continuum couplings. In this case, 
we found that the effect of the projectile bound excited state lpi/2 on the total, complete 
and incomplete fusion cross sections is quite weak (~ 10%). The couplings to the bound 
excited state lpi/2 only redistributes the complete fusion cross section (thick solid curve) 
between the elastic channel and this channel, the fusion contribution from the elastic channel 
being 1.7 — 3.4 times larger than from the bound excited state lpi/2 for the range of energies 
studied. We would like to note that these fusion excitation functions show similar trends 
as those obtained by Hagino et al. 0. We agree that complete fusion cross sections are 
strongly enhanced due to the couplings to the projectile excited states compared with the 
no-coupling case at energies below and just above the Coulomb barrier {Vb ~ 36 MeV for 
the elastic channel), whereas they are hindered at above barrier energies. 

In Fig. 2b, we show the effect of continuum- continuum couplings on the total and complete 
fusion cross sections of Fig. 2a. It is found that well above the Coulomb barrier, both total 
and complete fusion cross sections are suppressed compared with the no-coupling case, and 
enhanced well below the barrier. Just below the Coulomb barrier (34 MeV < Ec.m. < 
36 MeV), complete fusion cross sections are suppressed, but this is not the case for total 
fusion cross sections. In the present case, we found that couplings to the projectile bound 
excited state lpi/2 redistribute (dot-dashed curve) the complete and incomplete fusion cross 
sections, while the total fusion cross sections (dashed curve) remain nearly constant. With 
couplings to the bound excited state lpi/2, the contribution to complete fusion from the 
elastic channel is similar to the one from the bound excited state lpi/2 for energies below 
the Coulomb barrier, and 1.7 — 8 times smaller for energies above the Coulomb barrier. 

Fig. 3 shows incomplete fusion excitation functions (difference between the total and the 
complete fusion curves) for both cases presented above, namely in Figs. 2a and 2b. We can 
observe that continuum- continuum couplings significantly reduce the incomplete fusion cross 
sections (dashed curve). The case when the couplings to the bound excited state lpi/2 are 
not included is shown by the dot-dashed curve. 

From Fig. 2b and Fig. 3, it is observed that continuum-continuum couplings strongly affect 
the predicted total, complete and incomplete fusion cross sections. This implies that the 
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fusion dynamics strongly depends on continuum-continuum couplings. Since the short- 
ranged imaginary potential is well confined within the Coulomb barrier, we deduce that 
continuum-continuum couplings mainly reduce the fiux that penetrates the barrier, while 
couplings to the projectile bound excited state lpi/2 mainly redistribute, among the complete 
and incomplete fusion channels, the fiux that has already penetrated the Coulomb barrier. 

We have checked the convergence of reported fusion cross sections (total, complete and 
incomplete) with the size of the breakup subspace, and have found the following when 
couplings between all projectile excited states (bound-continuum and continuum-continuum) 
are included in the calculation: 

• The maximum energy of the continuum states (Fig.4): a maximum energy beyond 9 
MeV is needed to obtain converged results. With respect to the fusion cross sections 
for a maximum energy of 9 MeV (dashed curve), fusion cross sections for a maximum 
energy of 10 MeV (full squares) are changed by ~ 10% for energies around the Coulomb 
barrier. For energies well above the barrier, fusion cross sections are changed by ~ 
1.5%. 

• The density of the continuum discretisation (Fig.5): a density greater than 1.67 
bins/MeV is needed to obtain converged results. The same density is used for all 
partial waves in the continuum. With respect to fusion cross sections for a density 
of 1.67 bins/MeV (dotted curve), fusion cross sections for a density of 2 bins/MeV 
(dashed curve) are changed by ~ 6.5% for energies around the Coulomb barrier. For 
energies well above the barrier, fusion cross sections are changed by ~ 1.6%. 

• The number of partial waves in the continuum and potential multipoles (Fig.6): par- 
tial waves beyond /5/2, j-jji and potential multipoles beyond the octupole contribution 
(^=3) are needed to obtain converged results. With respect to fusion cross sections for 
continuum partial waves up to /-waves and potential multipoles K <?> (full circles), 
fusion cross sections for continuum partial waves up to (yf-waves and potential multi- 
poles K < (full triangles) are changed by ~ 8% (total fusion), ~ 100% (complete 
fusion) and ~ 3% (incomplete fusion), respectively, for energies around the Coulomb 
barrier. For energies well above the barrier, fusion cross sections are changed by ~ 
13%. The calculation including both continuum partial waves up to ^^-waves and po- 
tential multipoles 7^ < 4 (full triangles) is presently at the limit of our computational 
capability. 
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Fig. 7 shows experimental total fusion cross sections (full squares) for the similar system 
Be -|- ^"^Bi [jlO|, which should not differ too much from the reaction studied. By comparing 



converged total fusion cross sections for ^^Be + ^°^Pb (full stars), calculated within our 
model, with the experimental ones for ^^Be + ^°^Bi, it is observed that the converged total 
fusion excitation function does not reproduce the experimental one. They do agree with the 
experiment for energies around the Coulomb barrier, but underestimate the data by ~ 41% 
for energies well above the Coulomb barrier. 

A crude estimation of the effect of target excitations on the total fusion cross section 
has been done by (i) fitting the converged total fusion cross section in a single (elastic) 
channel calculation by finding an appropriate projectile-target real Wood-Saxon potential 
with an energy dependent depth and the geometry tq = 1.179 fm and a = 0.658 fm, and 
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then (ii) including the target excitations as in ref. [^. Such estimation reveals that the 
effect is quite weak. Fusion cross sections are increased ~ 1.28 times for energies around the 
Coulomb barrier, while they remain nearly constant for energies well above the Coulomb 
barrier. 

The experimental cross sections for ^^Be + ^°^Bi were obtained [TO as the sum of three 



channels: 5n+4n+fission. It was pointed out in ref. p8| that the 3n channel, expected to 
be relevant below the barrier, could not be measured and at the same time the fission cross 
section could have been overestimated. A new experiment is necessary in order to clarify 
the ^^Be fusion mechanism discussed in the present work. 

Summary and conclusions: Fusion cross sections calculated in the CDCC framework de- 
pend strongly on continuum-continuum couplings. We do not include transfer or inelastic 
channels of the target. Continuum-continuum couplings hinder total, complete and incom- 
plete fusion processes. Couplings to the projectile lpi/2 bound excited state redistribute 
the complete and incomplete fusion cross sections, but do not change the total fusion cross 
section. Results show that continuum-continuum couplings enhance the irreversibility of 
breakup and reduce the flux that penetrates the Coulomb barrier. A large breakup subspace 
is needed to obtain converged fusion cross sections. The converged total fusion excitation 
function does not reproduce the experimental one: converged total fusion cross sections 
agree with the experimental ones for energies around the Coulomb barrier, but underesti- 
mate those for energies well above the Coulomb barrier. A crude estimation of the effect of 
target excitations on the total fusion cross section reveals that it is quite weak. A new ex- 
periment seems to be necessary to clarify the ^^Be fusion mechanism discussed in the present 
work. The total fusion cross section is unambigously calculated in our formalism, but this 
is not the case for the complete fusion since the capture of all projectile fragments from 
breakup channels cannot always be distinguished from the capture of only one projectile 
fragment. 
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FIG. 1. Continuum discretisation used to define the energy bins. The central energies of bins 
are shown by full circles. 
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FIG. 2. Fusion cross sections as a function of the bombarding energy in the center of mass system 
for ^^Be + ^°^Pb. a) Include only couplings from and to the ^^Be bound states, b) Couphngs 
between all ^^Be excited states (continuum-continuum) axe included. See text for further details. 
The arrow indicates the Coulomb barrier for the elastic channel. 
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FIG. 3. Incomplete fusion excitation functions for both cases shown in Figs. 2a (solid curve) and 
2b (dashed curve). Incomplete fusion excitation function for the case of Fig. 2b, but couplings to 
the ^^Be bound excited state lpi/2 are not included (dot-dashed curve). 
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E (MeV) 

cm. ^ ' 

FIG. 4. Convergence of total (a), complete (b) and incomplete (c) fusion excitation functions 
for ^-"^Be + 208pj-, -^^j^j^ regard to the maximum energy of the -"^^Be continuum states included in the 
calculation. The s-, p- and d-wave continuum states for a density of the continuum discretisation 
of 2 bins/MeV, potential multipoles if < 2 and couplings between ^^Be excited states are included 
in the calculation. See text for further details. 
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FIG. 5. The same as in Fig. 4, but with regard to the density of the continuum discretisation. 
The maximum energy of the ^^Be continuum states is 9 MeV. 
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FIG. 6. The same as in Figs. 4 and 5, but with regard to the number of partial waves in the 
continuum and potential multipoles. The maximum energy of the ^^Be continuum states and the 
density of the continuum discretisation are 9 MeV and 2 bins/MeV, respectively. 
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FIG. 7. Converged total fusion cross sections for ^^Be + ^'^^Pb (full stars) are compared with the 
experimental ones [IC] for ^^Be + ^''^Bi (full squares). A maximum energy of the ^^Be continuum 
states of 10 MeV, continuum partial waves up to gr-waves for a density of the continuum discreti- 
sation of 2 bins/MeV, potential multipoles K < 4 and couplings between ^^Be excited states are 
included in the calculation. 
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